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Recent work by T. G. Philbin [l[ has provided a Lagrangian theory that establishes a general 
method for the canonical quantization of the electromagnetic field in any dispersive, lossy, linear 
dielectric. Working from this theory, we extend the Lagrangian description to reciprocal and non- 
reciprocal magnetoelectric (bi-anisotropic) media, showing that some versions of the constitutive 
relations are inconsistent with a real Lagrangian, and hence with quantization. This amounts to 
a restriction on the magnitude of the magnetoelectric coupling. Moreover from the point of view 
of quantization, moving media are shown to be fundamentally different from stationary magneto- 
electrics, despite the formal similarity in the constitutive relations. 
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I. INTRODUCTION 

In classical electromagnetism life is made much simpler 
by the introduction of e & fi. The microscopic current 
and charge densities can be forgotten, and in place of 
these we can deal with the macroscopic Maxwell equa- 
tions Indeed, for the purposes of further speeding 
up calculations, classical electromagnetism allows us to 
make artificial simplifications, such as the existence of 
media without dispersion, or loss. 

Quantum mechanics does not take to this description 
so easily. Although the quantization of the free electro- 
magnetic field can be found in standard textbooks Q, 
there has been some historical difficulty in quantizing the 
electromagnetic field within a dielectric medium. The 
field can be quantised in a fictional medium without dis- 
persion or loss [|[ , but as soon as dispersion is introduced, 
the procedure becomes awkward [5[. Moreover, it was 
not immediately obvious how an effective description of 
loss might be implemented in quantum mechanics. The 
classical field amplitudes ought to decay, whereas their 
operator counterparts must satisfy the canonical commu- 
tation relations uniformly throughout space, at all times. 

In fact these apparently distinct difficulties have their 
origins within a single physical effect. Dispersion and loss 
are the two sides of one phenomenon; the finite response 
time of a material to events that happened in the past. 
It is therefore through re-introducing some degrees of 
freedom associated with the medium that quantization 
may be carried out. 

Canonical quantization was achieved for a model La- 
grangian by Huttner and Barnett Q, who introduced a 
bath of harmonic oscillators to account for the dynam- 
ics of a uniform, dispersive and lossy dielectric. Sub- 
sequently, this model was extended to non-uniform di- 
electrics 0, 13 , and more recently, it was recognised that 
some model aspects of the theory could be removed, and 
that the theory could describe general features of non- 
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Most recently, a Lagrangian density was found in [l| 
that describes the electromagnetic field within any lin- 
ear magnctodielectric that satisfies the Kramers-Kronig 
relations. From this Lagrangian a Hamiltonian was de- 
rived that allowed a canonical quantization of the electro- 
magnetic field, and from this, formerly phenomenological 
results, such as the theory of Casimir forces (Lifshitz the- 
ory) have been given a canonical basis [l3[ . 

The approach here is to take seriously the Lagrangian 
in [l| . This Lagrangian represents an arbitrary linear ma- 
terial that automatically satisfies the Kramers-Kronig re- 
lations, as well as some general properties usually arising 
from thermodynamics (see section[H|). On top of this, the 
theory can be quantized, which suggests that this may be 
a more fundamental, and correct way to describe macro- 
scopic electromagnetism: the material degrees of freedom 
have returned, but only in the most minimal way. 

We examine extensions to the Lagrangian that describe 
the effects of magnetoelectric coupling, as well as time 
irreversibility (e.g. a medium in an external magnetic 
field). Here we understand magnetoelectric materials to 
include chiral media, moving media, and any other media 
where the constitutive relations arc of the form, D = 
e • E + x E b * an d H = fi -1 • B - x BE - E. Throughout 
what follows, a tilde over a vector denotes it being in the 
frequency domain. 

The aim is so that; (a) we may understand the origin of 
these various effects in terms of interaction terms within 
a Lagrangian; and (b) we may ask whether the require- 
ment of the existence of a corresponding Hamiltonian 
(from which we may quantize the field) places restrictions 
on the parameters within the constitutive relations. Wc 
find the answer to (b) is positive, a result which may clar- 
ify the apparent confusion over the restrictions placed on 
magnetoelectric media [4l[ (see [3413 )• This confusion 
between the various proposed restrictions is not negligi- 
ble, and is important for metamaterial design [lsL Il9|. as 
well as a possible route to a repulsive Casimir effect [2(| ■ 
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II. A SUMMARY OF THE LAGRANGIAN 
THEORY OF MACROSCOPIC 
ELECTROMAGNETISM 

We begin with a brief review of the basic features of 
the Lagrangian theory of macroscopic electromagnetism 
presented in [l[ . It is worth emphasizing that here we are 
concerned with the existence of a Hamiltonian that can 
be used to describe the interaction of electromagnetism 
with more general materials. We do not examine the 
subsequent quantization procedure, leaving this aspect 
for a future publication. 



1. The Lagrangian and the equations of motion 

The Lagrangian density of macroscopic electromag- 
netism can be motivated as follows. Firstly we have the 
familiar term associated with the electromagnetic field in 
vacuum, 



J2f p = | [E 2 - c 2 B 2 ] 



(1) 



where, E = - V0 - A, and B = VxA 

All media are dispersive, and hence via the Kramers- 
Kronig relations exhibit significant loss at some frequen- 
cies. When the medium does not depend explicitly on 
time, the Lagrangian must conserve energy, and this field 
energy lost from the dynamics of ([T|) must be transferred 
into another system. The response of the medium is as- 
sumed to have an arbitrary spatial dependence, so the 
additional system is proposed to be a reservoir of in- 
dependent harmonic oscillators that exists at all points 
in space, and at each point contains every possible fre- 
quency of oscillator, 
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X,., 



Yl-u 2 (xl + Yl) du> (2) 



In the Lagrangian density associated with the reservoir, 
d2J), two oscillators, X^ and Y u , are present at every 
frequency to account for the fact that we must distin- 
guish between a loss of field energy through the electric 
interaction (the imaginary part of e [42|), and through 
the magnetic interaction (the imaginary part of fi) . 

Finally, the field must be coupled to the bath of oscilla- 
tors in such a way that the classical macroscopic Maxwell 
equations give an extreme value for the action. It is found 
that such a coupling is given by, 



=^[NT E ' 
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/ a EE (cj)-X w dcj + B- / a± 
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(3) 

where the interaction with the material is determined 
via the coupling tensors, oc ee (uj) and a BB (u), which are 
assumed to be analytic functions of lu in the upper half 
plane, and vanish in vacuum. Note that all quantities 
appearing in (jT]{3]) are implicitly functions of position, 
and all of the fields arc also functions of time. 



The interaction, can be understood in terms of 
the local polarization, P, and magnetization, 7W, of the 
medium, 
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(4) 



With this notation, © takes the usual dipolar form, 
which would be expected from a local interaction with 
a neutral medium (see e.g. plj). 

As usual, the action is equal to a four dimensional inte- 
gral of the Lagrangian density, which in this case is given 
by Se =© + ©+©, 

S[<j>,A,X u ,Y u ] = J [C F +^ R +X,T}d 4 x (5) 

The remarkable features of © only become clear in the 
equations of motion, which are derived from finding an 
extremum of S. 

Using the usual field equations [22j], each of the oscil- 
lators is found to evolve according to, 



X u = -oj 2 X u + oc^ e (uj)-E 
Y L , = -uj 2 Y L , + al B (u)-B, 

with the field obeying the usual Maxwell equations, 

V-D = 

dD 



(G) 



V X H 



dt 



(7) 



where, D = e E + P, & H = B//j, - M, with P and 
M determined by the solution to ©, via (|1]). 

When the coupling, a EE , docs not depend on time, 
the simplest way to solve ^ is in the frequency domain, 
where, for example, X u = J X u (Cl) exp (— iVLt)dQ,/2ir . 
This leads to, 



a T E (uj).E(n) 



(w + + irj){u) — f2 — irj) 

+ 2tt [h xu 8(n -u) + h*J(fl + oj)} (8) 



where the h x (Y ) W are solutions to the homogeneous equa- 
tion. A similar relation also holds for Y u . We under- 
stand the transformation of (jSJ) into the time domain in 
the limit as, i] —> 0, where the choice of poles has been 
constrained by , X U (Q) = X w (— f2). The choice of the 
sign of rj corresponds to retarded versus advanced solu- 
tions. As we shall see, in this case the choice also amounts 
to the physics of absorption versus gain. 

Transforming © into the time domain, along with the 
corresponding expression for Y u , we insert these quan- 
tities into (0]), and find the general evolution of P and 
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M in terms of the field amplitudes and the h u . For the 
electric polarization, this gives, 



P ^ 



2^ 



du) 



oc ee (uj) ■ a^(w) • E(n)e 



-int 



±i 



dw 
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E 



( W ) 



( W 2 - ft 2 ) 



4w 



(9) 



where 'P' denotes the Cauchy principal value of the in- 
tegral, and the choice of sign comes from the choice of 
sign of r\ before the limit is taken. If we consider © in 
the case of sgn(?7) = +1, then it becomes clear that if we 
write, 

q ee ( bb)H • a EE(BB) (w) = ^3[x EE(BB) H] (10) 

with, Y B ,,„ B ,(ii)) interpreted as the electric (magnetic) 
susceptibility [43[ , then the Kramers-Kronig relations [2[ 
are automatically satisfied, 
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E (B B) \ 
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-dn, (ii) 



and that the electric polarization in @ can be written 
in the usual form, 



du 
2^ 



XeeM • E{u)e 



(12) 



where, Pq = a EE (uj) ■ h XLJ exp (-iuit)duj + c.c. is the 
'undriven' part of the polarization, that does not depend 
on the electromagnetic field. Similar relations hold for 
the magnetization, M, with the substitutions E — » B, 
and X — s* Y. It is striking that, with the identification 
given in (|10[) . the classical equations of motion arising 
from (O can represent electromagnetism within any ma- 
terial that may be characterised by linear electric and 
magnetic susceptibilities. 

The fact that the theory can be written in terms of a 
Lagrangian containing only first order time derivatives, 
from which we can derive a Hamiltonian, enables it to 
be quantized canonically; with Poisson or Dirac brack- 
ets becoming commutators. In quantizing the field, Pq 
and Mq are present by necessity: the classical ampli- 
tude, h u becomes an operator, related to the fluctuating 
'noise currents' already found to be necessary within the 
phenomenological theory (e.g. [H}). Furthermore, as the 
Kramers-Kronig relations, (fTTj) . arise from the dynam- 
ics of the reservoir of oscillators, it should be possible 
to apply this technique to quantize any theory of linear 
response. 



2. Restrictions on the constitutive relations, and the origin 
of loss versus gain 



If we define the coupling to the reservoir via ([10]) , then 
the susceptibility tensors are only consistent with the La- 



grangian if they are symmetric; a requirement usually 
arising from statistical physics Q, 24] , but automatically 
fulfilled here. Due to this symmetry of the susceptibility 
tensor there is also a well defined procedure for determin- 
ing a E E (B B ) from $s[x E E ( B b>] > although the result will not 
be unique [44|. The procedure also only produces a real 

O^EE (B B) 

if 9[Xee ( bb)] is positive definite. 
The lack of uniqueness of the coupling tensors plays no 
role in the classical theory, as a E E (B B) only appears lin- 
early in Pq, which may be given any form through a suit- 
able choice of h^. However, when the system is quantized 
the equivalent of the 'noise current' operator is given in 
terms of the operator versions of Po and Mq (with the 
h u now becoming operators with certain commutation 
relations) [l|, and therefore the choice of a EE(B B) will 
make a difference to the vacuum fluctuation of the polar- 
ization and magnetization of the medium. It is not clear 
whether this difference has any observable consequences, 
although we do not address this problem here. 

Another interesting feature is that the sign of the imag- 
inary part of X EE(BB) in (O is determined by a choice 
made in the dynamics of the reservoir. Therefore the 
distinction between loss and gain in this theory does not 
appear at the level of the Lagrangian, but in the bound- 
ary conditions imposed on the equations of motion of the 
reservoir. 



III. AN EXTENSION TO OTHER 
RESERVOIR-FIELD COUPLING TERMS 



As stated in the introduction, there are several reasons 
to investigate an extension of the Lagrangian theory of 
section [TT] to magnetoelectric media. There is a history 
of controversy regarding not only bounds on the mag- 



netoelectric susceptibilities, Xe 



but also the 



kinds of magnet oelectric coupling that are possible in 
principle [25l-l28| . Moreover, it is not clear to what ex- 
tent moving media and magnetoelectrics are equivalent, 
something which is important in the discussion of a 'fric- 
tional' component to the Casimir force [29l - [3lj |. Finding 
a Lagrangian that describes these materials and can be 
quantized may help clarify these issues. 

If there are macroscopic parameters that prove incon- 
sistent with a Lagrangian (or more precisely, a Hamil- 
tonian), the electromagnetic field cannot be quantized 
within such media in an obvious way, and we take this 
as a restriction on the material parameters. 

What additional coupling terms could be added into 
([21)? Lagrangians containing time derivatives of the fields 
higher than the first (which cannot be removed with a 
gauge transformation) do not have well defined canonical 
momenta, and consequently we cannot derive a Hamil- 
tonian. If we also consider the polarization and mag- 
netization to only depend upon the local values of the 
oscillators, X u & Y u , then we have the following gener- 
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TABLE I: Symmetry of coupling tensors under spatial inver- 
sion and time reversal. A plus sign indicates that the coupling 
tensor does not change under the inversion operation, and a 
minus sign indicates that the tensor is multiplied by —1. 



alization of ([4]), 



/3 BE . On the other hand, inversion symmetric, time irre- 
versible media, such as a dielectric in an external mag- 
netic field (a Faraday medium) must be described by E E 
& (3 BB . Finally, time-irreversible magnetoelectrics that 
also violate spatial inversion symmetry, such as moving 
media, or Tellegen media [32| must arise from a EB and 

We now examine the macroscopic theory that results 
from using the interaction Lagrangian with the polariza- 
tion and magnetization, ([TBI [T4"|) . 



A. The evolution of the oscillator amplitudes 

With the interaction Lagrangian defined by (|13J414p as 
in ((3]), the equations of motion for each of the oscillators 
in the reservoir arc now, 



and, 



£*EE + /3 EE r^ ) • X u 



Q EB+ (3 E B J m Ytt 



du (13) 



and, 



a T - 

E E 



d 



:E dt 

T 



• E 

-B T - 

U BE rBE^ 



B (15) 



M 



O-B B + /3 B Bgj ] ' Y ^ 



"be +^be^ I • X U 



dbJ (14) 



Table U shows the properties of the coupling tensors 
that arise from assuming that the value of the action 
is separately invariant under an active time reversal 
(t — > —t) and space inversion (x —> —x) of the fields; 
e.g. Xj{x,t) -> X u (x,-t) k X u (x,t) -> X u (-x,t). 
In obtaining table HI we have assumed that the X u & Y u 
oscillator amplitudes obey the same symmetry properties 
as the electric and magnetic fields, respectively. 

If a coupling tensor is non-zero and changes sign un- 
der time reversal, then the medium possesses an intrin- 
sic time irreversibility (e.g. due to an external magnetic 
field, or motion). Meanwhile a change of sign under spa- 
tial inversion shows that the medium possesses a certain 
handedness (e.g. chiral media). 

Table U therefore shows that magnetoelectrics violat- 
ing spatial inversion symmetry, while exhibiting time- 
reversibility, such as those constructed from chiral inclu- 
sions, must be described by the coupling terms /3 EB and 



-u 2 Y u 



/3bb^t -B 



. T 



x 9 



■ 3 BB -).E (16) 



In Fourier space, we therefore obtain, 



(u> — Q — irf) (lj + Q + irf) 



+ mal^j -b + (c 



(lu — — irf) (cj + £1 + irf) 



E 



(17) 



with the poles dealt with as in section HH and the ho- 
mogeneous parts of the solution omitted (c.f. ©). Note 
that from now on we consider the case of media with loss 
rather than gain. 

Transforming (|17p into the time domain gives us the 
final expressions for the evolution of the reservoir. The 
Xu oscillator obeys, 
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"L+'^EB - E + k B +^BB B 



e- laJi + ^e- lwt +c.c., (18) 



and the 



"bb + ^bb ■■»+ aJ B +in/3 EB 



(a; 2 - fi 2 ) 



4w 



c% B +iu0£ B )-B+(al B +iLj0T B )-E e~ iut + h^e** + c.c. (19) 



Equations (jT51 IT^j) now determine the evolution of P 
M in terms of the held amplitudes. 



B. The polarization and magnetization in terms of 
the held amplitudes 

Inserting (fTSITT^ji into (p^TH)) gives us a polarization 
and magnetization of the form, 

XeeH • -E + XebM • B\ t~" 



where, 



p _ p 

/OO 7 
X BB H-B + X BE (^)-£;] e"^ — 

(20) 

where the electric-electric and magnetic-magnetic sus- 
ceptibilities are given by, 

, A EE • A EE (^ ,Q)+A EB -A EB (q;,0) 
(fi 2 -w 2 ) 

+ A EB • A 

A BB • A BB (w,r>) + A BE • A BE (w,r!) 



x EE (w) =p / <m 

10 

in 



^[ a ee • A EE (w,w) + A EB • A EB (w,w)] 



/•OO 

x BB M =p / dn 

Jo 



(fl 2 -UJ 2 ) 



+ ^j[ A bb ■ A BB (w,w) + A BE • A BE (w,w)] 

(21) 

and the magnetoelectric susceptibilities are, 

A EE • A BE (u; ,r2) + A EB ■ \j B (uj,n) 

(n 2 -cj 2 ) 

+ A E B • A. 

^bb ' A EB (w,ft) + A BE • A EE (w,fi) 



Xbb (w) =P / dfl 
Jo 

ITT 



— [A BE • A be (w,cj) + A EB • A BB (w,w)] 



p OO 

x BE M =p / dn 

Jo 



(n 2 -Lu 2 ) 



— [A BB • A EB (u;,w) + A BE • A BB (w,w)] 



(23) 



(22) 



A EE (u;,Q) = a EE (0) - iw/3 EE (fi) 
A BB (w,0) = a BB (^) - iw/3 BB (fi) 
A eb (cj,0) = a EB (^) - iw/3 EB (f2) 
A BE (w,0) = c*be(^) - iw/3 BE (Q) 



The physical interpretation of the dependence of (|23|) on 
two frequencies, w and f2, is that in general the medium 
responds to the field amplitudes and the rate of change 
of the field amplitudes. Therefore, in the integrals within 
(|2~T1l2"2"l) we have factors of u> that arise from a linear 
response to a change in the field (e.g. a polarization 
due to the EMF from a changing magnetic field), and 
factors of f2 that relate to the dispersion due to the finite 
response time of the medium. 

With the results ([2~T1 122]) . we can see that the suscepti- 
bilities arc divided into parts related to dissipation, and 
lossless response, the connection between the two being 
an interesting generalization of the Kramers-Kronig re- 
lation (fTTj) . 

For the electric-electric and magnetic-magnetic sus- 
ceptibilities, the dissipation is determined by the anti- 
Hermitian part of the susceptibility (c.f. 0), while the 
magnetoelectric dissipation is determined by x BB -^L 
(as can be verified through substituting the magnetoelec- 
tric constitutive relations into §80 of [3|). It is interesting 
that these two susceptibilities are automatically related 
in such a way that a non-zero % BB necessitates a non- 
zero Xbe- 

The undriven parts of the polarization and magnetiza- 
tion in (|20p involve a coupling between the two reservoirs, 



P OO 

5 o = / [A EE (w,u;) • + X EB (u,u) ■ h YU ]e~ %UJt du 
Jo 

+ c.c. (24) 
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and, 

/>oo 

M = [\ BB (uj,uj)-h XUJ + X BE (uj,uj)-h VUJ ]e~ tult duj 
Jo 

+ c.c. (25) 

Having now established a general form for P & M, we 
develop the physical interpretation of the various cou- 
pling tensors within the Lagrangian. 



C. Physical interpretation 

The constitutive relations (|2T)|) are, as anticipated, 
those of a magnetoelectric. However, this is a very gen- 
eral form, where the material neither obeys time reversal, 
nor spatial inversion symmetry. We now consider three 
special cases of ([2T][22]); time reversible, odd parity me- 
dia; time irreversible, even parity media; and media that 
exhibit a static magnetoelectric response. 



1. Time reversible, odd parity media 

If we demand that the medium be time-reversible, 
then, from table, HI a EB = a BE = (3 EE = /3 BB =0, and 
we obtain the frequency domain constitutive relations, 

D = e • E + iujK ■ B 

H = H 1 • B + iuiK T • E 

where, e = e Q 1 3 + x EE , A 4 " 1 = A^ 3 ^ - X BB , and 
2wQ ! (k) = a EE • /3g E — /3 EB • c* BB , and the real and imag- 
inary parts of k are connected by (fTTjl . Consistent with 
the time reversibility of this situation, the permeability 
and permittivity are symmetric tensors. 

These are the constitutive relations of an anisotropic 
chiral medium in Boys-Post form [33| . Notice the neces- 
sity of having n T appearing in H , versus n appearing in 
D. No other choices of k appear to be consistent with 
the theory. The prefactor of lu in the magnetoelectric 
coupling is crucial, as it has the consequence that in the 
static limit (lu — > 0), unless k diverges, the medium no 
longer has a magnetoelectric response. 



2. Time irreversible, even parity media 

Applying table U in this case means that ct EB = a BE = 
f3 EB = f3 BE = 0, and the magnetoelectric susceptibilities, 
(|22"jl vanish. We are then left with, 

D = eE 
H = H 1 B (26) 

where e = e l 3 + X EE , and fi 1 = ^3 - X BB - 



The consequence of having broken the time reversibil- 
ity of the medium is that at frequencies where the loss 
is negligible, e = and fi = fv . The fact that these 
tensors become Hcrmitian is consistent with the gener- 
alized principle of the symmetry of kinetic coefficients in 
the case when time reversal symmetry is broken, and the 
medium is without loss [U, [24| . Again, results from sta- 
tistical physics emerge from this Lagrangian description. 
As lu — > 0, e and /x become symmetric tensors. 

One physical example of a medium which would be 
described by a /3 EE term would be an ordinary dielectric 
in an external magnetic field. For frequencies where such 
a medium is without loss, e is a Hermitian tensor 0]. 



3. Media exhibiting a static magnetoelectric response 

In the limit when the electromagnetic field becomes 
static (lu — > 0), all of the quantities in (f2"5| become 
real, and only the a coupling tensors play a role. The 
reason for considering this limit is that here we isolate 
the terms in the Lagrangian that describe the physics 
of media in motion, and distinguish Tcllcgcn from chi- 
ral media (something that has caused controversy in the 
past [2(1 [2|| ) . The static field constitutive relations are, 

D(lu^0) = €-E + Xeb -B 

H(uj -+ 0) = /x- 1 • B - xl B ■ E (27) 

Observe that, if a medium is to have a static magneto- 
electric response, then this must be characterised with 
the coupling terms containing the a EB and ct BE tensors. 
In this limit, e and \x become symmetric tensors, and 
the magneto electric coupling, x EB must appear with a 
minus sign and a transpose in H versus D. It is incon- 
sistent with the Lagrangian to suppose that the static 
magnetoelectric coupling can have any other form. 

The constitutive relations, ((2"7]) . when extended to 
arbitrary cj, are consistent with the kinds of magneto- 
electric coupling required for both moving media and 
Tellegen media jU, HH, HI[ • From the point of view of the 
matter-field coupling terms in the Lagrangian descrip- 
tion, there does not appear to be any contradiction in 
assuming that through violating both parity and time 
reversal symmetry, such media could be constructed in 
the laboratory. 

This may well be true for Tellegen media, however, for 
real moving media there is a subtlety that means that 
the magnetoelectric coupling in (|2"T|) docs not contain 
all of the physics of electromagnetism interacting with 
a medium in motion. As shall be shown in section IIV1 
the reservoir Lagrangian density, Jz? R also has to be al- 
tered in this case, and it is not obvious how one would 
engineer a medium where the loss mechanism works in 
such a peculiar way. 
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D. Further restrictions on the constitutive relations 

From the above discussion, it appears that the general 
susceptibility tensors, ([2"Tll2"2")) encompass all known mag- 
netoelectric constitutive relations, and naturally restrict 
the relationship between % EB and Xbe- 0\rr point of 
view is that this is the correct description, as the electro- 
magnetic field may be quantised within such a formalism. 
This point of view is bolstered by the fact that some of 
the expected restrictions from thermodynamics have also 
arisen along the way. 

This formalism puts further restrictions on the suscep- 
tibilities, due to the fact that (j2"Tll2"2"j) are related to one 
another. For example X EE & X BB together contain the 
same eight coupling tensors as x EB . Therefore it is not 
possible to choose the magnetoelectric coupling in a way 
that is independent of the value of the electric-electric 
and magnetic-magnetic susceptibilities. We now proceed 
to work out the implications of this relationship. 

The dissipative part of the susceptibilities can be sum- 
marized as follows, 



L-^EE XeeJ 



rt EE E E ' A EB E B 



Xb b] 



ITT 

~J~ [^eb Xbe] 
ITT 



-^B B * A B B + A B E * A B E 



Ae e * A B E -f- A E b * A B B 



(28) 

The problem of finding the constraints on the compo- 
nents of the magnetoelectric susceptibilities versus the 
electric-electric and magnetic-magnetic ones is now one 
of linear algebra. There are four complex matrices, and 
we have to work out how the components resulting from 
multiplying them together in one way are related to mul- 
tiplying them together in another way. In appendix 
this is calculated, and the following restriction on the 
susceptibility tensors is obtained, 



* f 1 

A.B El 



< 



[Xe E Xe e] 5 



[Xb b Xb b] - 



(29) 

If the Lagrangian with the coupling (fT3"l[L"f| is taken 
to describe the most general kind of magnetoelectric 
medium, then every such medium should exhibit dissi- 
pation that satisfies (f2T))) in order to be consistent with 
quantization. 

It is worth examining (|29p in a few specific cases. In 
the case of chiral media, using the notation of section 
MI C II the inequality becomes, 



, 2 C> 



which in the isotropic case is identical to the result of [16[ 
(0). 

For the case of magnetoelectrics in the uj — > limit, 
the inequality is neither that of [l4[ nor [l5| , 

^[XebLj — ^[Aee]m^[Xbb]«' 



The restriction (|2T?|) affects the dissipative part of the 
susceptibilities at all frequencies. Of course, this restricts 
the non-dissipative parts in some way as well, via the 
Kramcrs-Kronig relations. However, for a given fixed 
frequency, the Lagrangian does not appear to place re- 
strictions on the non-dissipative parts of (|2"T1 12"2"]) . This 
is contrary to the inequalities that are often used in the 
literature, which do not seem to fully treat dispersion or 
loss. 



IV. THE LAGRANGIAN FOR A MOVING 
MEDIUM 



It was established in section MI C 31 that the magne- 



toelectric coupling terms involving a E 



within the 



Lagrangian can in principle reproduce the constitutive 
relations for a moving medium (with the f3 tensors equal 
to zero). This section is motivated by a recent discus- 
sion regarding the existence of a frictional component to 
the Casimir force [3(1 HH , where it has been pointed out 
that if a magnetoelectric could perfectly mimic a moving 
medium, then the frictional force could extract work from 
the vacuum. There is either no frictional force, or a mag- 
netoelectric cannot perfectly reproduce the physics of a 
moving medium. Here we argue that, from the point of 
view of a Lagrangian description, a magnetoelectric cou- 
pling alone is not sufficient to mimic a moving medium. 

Suppose that in the comoving (primed) frame we have 
a medium that can be described via e and fj, alone, as 
was assumed in section ILTI 



B'- / a4 B ( w '). (30) 
Jo 



Now consider the lab (unprimed) frame, where the 
medium is in uniform motion. Without loss of generality, 
we can assume the motion is along the x axis, V = V x x. 
The form of the Lagrangian density associated with the 
free field is unchanged in terms of the field strengths, as 
it is a scalar formed from F^F^ . However, Jr? INT and 
J2? R will not take the same form in terms of the oscilla- 
tor amplitudes in both frames. To find the form of the 
Lagrangian density in terms of lab frame quantities, we 
begin by transforming the field strengths, 

E x = E x B x = B x 

E' y = 1 {E y -V x B z ) B' y = 1 {B y + V x E z /c 2 ) 
E' z=1 (E z + V x B y ) B' z=1 (B z - V x E y /c 2 ) (31) 

with, 7 = (1 - V 2 /c 2 ) -1 / 2 . Inserting these into ([3D]) 
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yields an interaction Lagrangian, 



p oo 

Jz? INT = E • I [ol ee {lo) ■ X u + q eb (w) • YJ[ du 
Jo 

POO 

+ B- [a EB ( W ).y u |a BE ( W ).IJi (32) 
Jo 

with the following coupling tensors, 

a EE {uj) = A • a' EE (uj) a EB (uj) =7^X c^bM/c 2 
a BB (w) = A • a BB (ui) a BE (uj) = -jV X a EE (uj) 

(33) 

where A = diag(l,7, 7). When the medium is non- 
uniform in the rest frame, the coupling tensors become 
functions of time in the lab frame (e.g. a EE (x') = 
a EE (7(a; — V x t))). Thus far the coupling tensors have 
been assumed to be independent of time. 

The primes on the oscillator amplitudes and frequency 
have been dropped in obtaining (|3"2"|) from (f5D|) . This is 
because, as explained in the introduction, these oscillator 
amplitudes are an unobservablc accounting device for the 
lost field energy. We only have to make sure that their 
dynamics and coupling to the field are properly described 
in terms of lab frame coordinates and fields. The ampli- 
tudes themselves cannot be observed in either frame, and 
so there is no useful meaning in transforming them. 

As initially anticipated, in terms of the interaction La- 
grangian, a moving medium falls into the category of a 
time-irreversible medium without spatial inversion sym- 
metry (see table |T|. In both types of inversion the velocity 
changes sign, which represents different medium. Notice 
that the transformed interaction Lagrangian, Q32[) auto- 
matically contains the Aharonov-Casher interaction (3(| , 
V • M X E/c 2 , which has a subtle origin in the multipolar 
expansion [2l| . 

Despite the familiar interaction Lagrangian, moving 
media are fundamentally distinct from stationary mag- 
netoelectrics. This is due to the behaviour of the reservoir 
part of the Lagrangian density. In the rest frame this is as 
in ([2]), but with primed quantities. However, the deriva- 
tives of the oscillator amplitudes with respect to time 
in the rest frame, dX'^jdt' & dY'^t/dt', are not equal 
to the derivatives with respect to time in the lab frame. 
To describe the dynamics of the reservoir correctly in 
the lab frame, the time derivative must be transformed; 
d/dt' = "f(d/dt + V • V), and the reservoir part of the 
Lagrangian density becomes, 



7 2 (^ + (V-V)X U 



2 Jo V dt 



fdY u 



V dt 



(V-V)Y UJ ) -u 2 (Xj + YJ) }dw 



(34) 



where the primes are again dropped from the amplitudes 
and frequency, for the reason described earlier. This is 



not the same modification to the reservoir Lagrangian 
that was made in [l2| , although as shown below our form 
does produce the correct constitutive relations. 

We propose that the sum of (P), & flM]) repre- 

sents the Lagrangian for the description of a medium in 
uniform motion, when, in the rest frame this medium 
can be described by the tensors e and fi [46| . An exten- 
sion to non-uniform motion is possible through consid- 
ering a local rest frame Lagrangian at each point in the 
medium, however this case is not considered here. No- 
tice that the key feature of (fM]) is that the reservoir is 
fundamentally altered, even in the absence of the field. 
This is encoded within a coupling between neighbouring 
oscillators through the terms, (V • V)X W & (V • V)>^. 

To show that Jt? = dU) + d32) + (GEO) is the correct 
Lagrangian, we examine the polarization and magneti- 
zation of the medium that arises from the equations of 
motion. Examining (|32[) shows that the polarization and 
magnetization are given by, 

POO 

P= [a EE (w, x,t) ■ X u + at EB (uj,x,t) • Y u ] duj 
Jo 

POO 

/ [a BB (uj,x,t) ■ Y u + a BE ((j, x,t) ■ X u ] du 
Jo 

(35) 



M 



To find the quantities in ([33)1 in terms of the fields, we 
solve the equations of motion of the oscillators, which are 
now, 

^ (Ft + V ' ^ = + • E + «BE • B 



1 (dt + v ' V J YiJ= YbJ + BB ' B + EB ' E 

(36) 

The dynamics of (f3T>|) clearly demonstrates some kind 
of spatial dispersion. Assuming a uniform medium, we 
write p6p in Fourier space to find X^, 



X u (k,Sl) 



"beM 



B(k,n) 



(oj - Q' - irj) {ui + fl' + irj) 



2-7T7 



h^{k)$(u-n') + Kj-k)s(Lu + n') . (37) 



where SI' = — V • k), sgn(rj) = +1 (a moving 
medium with loss), and the result for Y ^ is obtained 
from interchanging the subscripts and fields, E B, & 
X <r> F.Note that spatial dispersion occurs in the de- 
nominator, which changes the positions of the poles. In 
short: the modified dynamics of the reservoir represents 
the physics of the Dopplcr effect, which is a very special 
kind of spatial dispersion. We emphasise that for this 
reason, this medium is not entirely equivalent to a sta- 
tionary one with a magnetoelectric coupling such as that 
of section IIII C 31 When (JHT)) is inserted in along 
with the corresponding expression for Y u , we obtain the 
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following polarization and magnetization vectors, 

r d 3 k r°° r 
P = P + — 3 / {[x BB (u,')-E(k,«,) 
Jr3 (2ir) Jv-k 1 



X EB (u')-B(k,u,j\i 



^{h-x — uit) 



z.c.ldw (38) 



and, 



d 3 k 



V-k 



{[ X ^')-B{k,w) 



/r3 (2tt) j 

+ Xbe ("') ■ E(k,u;)]e^ k - X -^ + c.c.}dw (39) 

where uj' = ^{ui — V-k), and the real and imaginary parts 
of the susceptibilities are related by (fTT|) . with uj — > u>'. 



The first thing to notice about ([38|) & (|39|) is that in 
the lab frame, some of the positive rest frame frequen- 
cies (uj' > 0) appear as negative frequencies (uj < 0). 
The response of the medium to a constant field, uj' = 0, 
also appears at a finite frequency in the lab, uj = V • k. 
These peculiar features are a consequence of the Doppler 
effect, which has arisen from the modified reservoir dy- 
namics encoded in (|34p . and would not have occurred 
for any stationary magnetoelectric. The undriven part 
of the polarization and magnetization also exhibits this 
behaviour, for example, 

f d 3 k f°° 
P o^l 7T~^3 / [a EE (uj') -h XUJ ,(k) 
y R3 (2n) 6 J v . k 

+ a EB ( W ')-V(fc)]e' M ^ + c.c. (40) 

The susceptibilities are given in terms of the coupling 
tensors as in section [III C 31 Inserting (|3"3"|) into (|2"Tl - |2"2")l . 
and identifying the rest frame susceptibilities according 
to (|10p . we find the following transformation formulae, 



Xe b 
Xb b 



— v x x'bb- a 



A-x' BB -A-j 2 V x X ' EE xV 



Xb e T 



V X Xee • A + -jA-Xbb x v 



which reduce to the well known first order in V j c results 
when the medium is isotropic (x' EE = Mx'ee & Xbb = 

D(k,uj) = e'(u;')E(k,u>) + (n'V) - i) V X H ^ w ^ 

<r 

B(k,u) = u,'(u')H(k,u) - (n'V) - y V * E ( k > u ) 

c 2 

(41) 



where, e'(w') = e + x' EB (o/), m' V) = /i 1 - Xbb( w 0; 
andn' (w')/c 2 = e' (uj')h' (u)'). Therefore, the La grangian 
density, _Sf = (P) + © + © reproduces the correct 
macroscopic Maxwell equations for a medium in motion 
with any e and fi that satisfy the Kramers-Kronig rela- 
tions. 

From the point of view of the quantization of macro- 
scopic QED in terms of a ficticous bath of oscillators, it 
may seem obvious that a moving medium should not be 
equivalent to a stationary magnetoelectric: a harmonic 
oscillator is not a rclativistically invariant system. How- 
ever, if we do not utilise the bath of oscillators in the 
quantization, then it is not obvious that the two are in- 
equivalent. In regimes where dispersion is negligible, the 
constitutive equations take the same form, and this ap- 
parent equivalence has led to a paradox in the theory 
of the Casimir effect [3(| , which we have hopefully shed 
some light on. 



V. CONCLUSIONS 

We have shown that a natural generalization of the La- 
grangian in [l[ provides a general description of magneto- 
electrics, and media exhibiting Hermitian e & fj, tensors. 
It has also been established that this generalization is 
only consistent with quantization (i.e. the existence of a 
Hamiltonian) with materials where the magnetoelectric 
coupling satisfies (|29[) . which reproduces some inequali- 
ties previously derived from different physical arguments. 
We propose to be an accurate restriction on the mag- 
nitude of the magnetoelectric coupling. 

It has also been shown that as far as the Lagrangian de- 
scription is concerned, moving media are not equivalent 
to stationary magnctoclectrics. The coupling between 
the field and the reservoir can indeed be the same in the 
two cases, however, the reservoir must be represented by 
([3~4"| rather than ^ to fully account for the physics of 
the Doppler effect. 
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Appendix A: A derivation of the magnetoelectric 
inequality 



Suppose we have four 3x3 matrices, a, (3, 7, & 6, all 
containing complex entries. We write these matrices in a 
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compressed notation as follows, 
an "12 "13 

tt21 "22 "23 
^"31 "32 "33 j 

where, <X; = (an, a.i2, aa). We now seek to find the 
relationship between the following matrix products, 



A 
v 
a 



a • a' 
a-S" 1 



-7-7' 
S ■ S^ 

7-/3 f 



Expanding out these matrix products, we find that their 
elements can be written as follows, 

Ay = o, • a* + -y l ■ 7* 
Uij = Pi ■ p* + Si • S* 
ciij <», • S] 7 ; • p* 

We now examine the matrix element, er.y and look to 
write it in terms of the diagonal terms, A„ & Vjj. Multi- 



plying together both these diagonal terms, we obtain, 



An Vjj 



|o*|W + l«<lW + l7ilW 



17/l^f 



(Al) 

We compare this with the absolute square of the element, 



'13 , 



|cr,;,'| 2 = la, • S*\ 2 



\ Tr p*\ 2 + 2dt [( ai .S*)( Tr p*)] 



A2) 

Applying the Cauchy-Schwarz inequality [4(|, |a;| 2 |y| 2 > 
\x-y\ 2 , to (O 

we thus obtain, 



< \OLi 



or, from (| Al I) . 



|7 i | 2 |/3 i | 2 +2|a i ||^||7 i ||/3 j 



So finally we find the following inequality must be satis- 
fied by the matrix elements, 



T ij I fl Ajj Vjj . 



